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Abstract 

We systematically investigate the possibility of finding CP/T violation in the r 
sector at Tau-Charm Factory. CP/T violation may occur at r pair production process, 
expressed as electric dipole moment, and at tau decay processes. By assuming that 
electric dipole moment as large as 10 -19 e-cm and CP/T violation effect orignating 
from r decay as large as 10 -3 are observable at Tau-Charm Factory, we studied all 
the possible extensions of the SM which are relevent for generating CP/T violation 
in t sector. And we pointed there are a few kind of models, which are hopeful for 
generating such CP/T violation. For these models we consider all the theoretical and 
current experimental constraints and find that there exists some parameter space which 
will result in a measurable CP/T violation. Therefore we conclude that Tau-Charm 
Factory is a hopeful place to discover CP/T violation in r sector. 
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1 Introduction 



The origin of CP violation has remained an unsolved problem since the discovery of CP 
violation in K meson system a quarter ago[|IJ. Although the observed CP violation in K 
meson system can be accommodated in the standard model (SM) of electroweak interactions 
by virtue of a physical complex phase in the three by three Cabibbo-Kobayashi-Maskawa 
matrix (CKM) ||, it is not clear if CKM mechanism is really correct or the only source 
for CP/T violation ||. To verify CKM mechanism one needs not only the information 
on K meson mixing and decay but also that from the B meson system or other systems. 
The main physical purpose of B factory is to test the CKM mechanism. However even if 
CKM is the correct mechanism to describe the CP violation in K and B meson mixing and 
decay, it is not necessary that the CKM matrix is the only source of CP/T violation in the 
nature ||]. As pointed out by Weinberg 0, unless the Higgs sector is extremely simple, it 
would be unnatural for Higgs-boson exchange not to contribute to CP/T non- conservation. 
CKM matrix may explain the observed CP violation in K meson system and possibly the 
CP violation in B meson system, while other new sources of CP/T violation may occur 
everywhere it can. In fact there are some physical motivations for people to seek the new 
sources of CP/T violation. One motivation is from strong CP problem in the SM ||. For 
most of the scenarios to solve this problem they need more complex vacuum structure and 
therefore new CP non-conservation origin. Another motivation is from cosmology, most 
astrophysical investigation shows that the additional sources of CP violation are needed to 
account for the baryon asymmetry of universe at present 0. The third motivation is from 
super symmetry. Even in the minimal supersymmetrical standard model (MSSM), there 
are some additional CP non-conservation sources beyond the CKM matrix ||. Now the 
question is at what places the possible new CP/T violation effects may show up and what 
is the potential to search for those effects. In this work we are going to study systematically 
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on the possibility to find new CP/T violation effects at Tau-Charm factory (TCF). 

TCF is a very good place to test the SM and search new physics phenomena because of its 
high luminosity and precision ||. Especially the r sector is a good place to seek for non-SM 
CP/T violation effects because in the SM CP violation in lepton sector occurs only at multi 
loop level and is way below any measurable level in high energy experiments, only non-SM 
sources of CP/T non-conservation may contribute and another reason is that r has abundant 
decay channels with sizable branching ratio, which can be used to measure CP/T violation. 
Furthermore, the production-decay sequences of r pair by electron-positron annihilation 
is also favored. The reason is as the following: (i) r pair production by electron-positron 
annihilation is a purely electroweak process and can be perturbatively calculated; (ii) For the 
unpolarized electron-positron collision, its initial state is CP invariant in the cm. frame; (iii) 
when the electron and/or positron beams are longitudinally polarized, the initial state is still 
effectively CP even, which presents extra chances to detect possible CP violation. To detect 
the possible CP/T violation, one can either compare certain decay properties of r~ with 
corresponding CP/T conjugations, or measure some CP/T-odd correlation of momentum or 
spin of the final state particles from r pair decay. These CP/T violating observables can and 
should be constructed model independently, since normally in non-SM these observables are 
not well predicted due to the complexity and many free parameters. The sensitivity of the 
experimental measurement on the possible CP/T violation is determined by the sensitivity of 
the measurement on momentum, spin or other physical quantities of the final state particles, 
from them the physical CP/T violating observables are constructed. The better one can 
measure these quantities, the momenta, for example, the smaller the CP violation phase 
can be reached. In TCF, it is to expect about 10 7 r pair in one year, and the precision of 
measurement on kinematic parameters at 10~ 3 . The statistical and systematic error can be 
around or below this level. Therefore generally a CP/T violation phase as small as order 
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of 10~ 3 can be reached at TCF |J. In a non-SM the CP/T violating phase may appear in 
various stages of the process of production-decay chain, e + e~ — > t + t~ — > final particles. 
We sort them in three cases; (i) CP/T violation is generated in the tree level production 
process, e + e~ — > ^,Z,X — > t + t~, where X is some new Higgs or gauge bosons, CP/T 
violating phase appears either in the propagator of X or in the coupling to lepton pairs, 
and the simplest possibility is X being a neutral Higgs in two or multi-Higgs doublets 
model. In this case the size of CP/T violation is proportional to the interference between 
the X exchange and 7, Z exchange processes. Unfortunately for X being Higgs doublet 
this interference term is proportional to the initial and final states fermion masses m e m T 
as a result of chirality conservation. This factor along contributes a suppression factor 
m e /m T ~ 3 x to all CP/T violating observables in this kind of processes at TCF 
besides other possible suppression factor, like the large mass of X, small coupling between 
X and leptons. We conclude that it is hopeless to search for CP non-conservation from the 
tree-level production process at TCF. (ii) CP/T violation is also generated at production 
stage, but through loop level. The most hopeful cases are that there may exist large electric 
or weak dipole moment (EDM or WDM) for r lepton, i.e. there are sizable CP/T violation 
phase at the vertex r~ — 7, Z — t + . For this situation the new physical particles beyond SM 
only appear as virtual particles through loops and the size of CP violation is proportional to 
EDM or WDM and is not suppressed by other factors, so the point is just whether EDM or 
WDM of r is large enough to be observed. Generally the Lagrangian describing the CP/T 
violation in r pair production related to EDM and WDM is 

L CP = -l/2tfa^ l5 r[df(q 2 )F, u + d^{q 2 )Z, u ] (1) 

F^ v and Z^ are the electromagnetic and weak field tensors. The momentum transfer at 
TCF is around 4 GeV, and in LEP experiments it is around the mass of Z boson. Therefore 
at TCF we expect the contribution from WDM is a factor of ^pp — 2 x 1CT 3 smaller than 
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the contribution from EDM, if EDM and WDM at the same order of the magnitude. On the 
other hand the EDM term is less important at LEP energy. That is the reason why the LEP 



data constrain more strictly on WDM than EDM of r |T(], We will neglect the WDM 
contribution from now on in this work, (iii) It is possible that the CP/T violation phase is 
small in the production process but it is relatively large in the r pair decay processes. The 
processes like r to neutrino plus light leptons or hadrons through some new bosons exchange 
at tree level can contribute significantly to CP/T violation observables. Obviously in this 
situation any CP/T violation effect from loop level is negligible, since any loop effect is 
at least suppressed by a factor j^jfi, where M is the mass of some new heavy particles 
appearing in loops. This factor is smaller than 10~ 4 if M is heavier than about 20 GeV. 

Now let us recall that how one detects CP violation in K meson decays: One measures 
the partial widths for a decay channel and compares it with that for the corresponding 
CP-conjugate decay process. Underlying such a philosophy is the interference between a 
CP violating phase and a CP conserving strong interaction phase, i.e. CP violation effect 
is only manifested in the process with strong final state interaction. To observe possible 
non-CKM CP violation effects in tau decays, however, one has to invoke new methodology 
in the most cases. The basic reason is that both in production vertex of r pair (EDM of 
r) and in some tau decay channels (like pure leptonic decay, m/, pv decay channels etc. ), 
there is no strong interaction phase, caused by hadronic final state interaction, to interfere 
with possible CP violating phase. So far some efforts have been made to investigate the 
CP/T violation effects in TCF. Mainly those work are trying to find various ways to measure 
possible CP/T violation. The simple and very useful method is to construct observables 
which are CP/T-odd operators being made from momenta of final state particles coming 
from r pair decay or polarization vector of the initial electron (or both electron and positron) 
beam | T2]j . These operators can be used very conveniently to test any CP/T violation from 
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either EDM of r lepton or from the decay of the r pair without much model dependence. 
Some of the operators are constructed by considering the reactions 

e + (p) + e -(-p) — r+ + r~ -> + + X (2) 

in the laboratory system, where A(B) can be identified as a charged particle coming from 
the t - (t + ) decay Some CP/T-odd operators (so CPT even, we will not consider CPT-odd 
operator in this work since it is certainly much smaller violation effect) can be expressed as 
following [T2[ 



Ot=p 



TV = {q + -q-y.!fe£Z + (i~j), (3) 



\q+X 0_|> 

(q+xq-) j 

\q+xq-\ 

where p,q denote the unit momenta. If the initial electron and/or positron beams are 
polarized, one can construct some more observables making use of the initial polarization 
vector. For example a T violating operator 

_ q+ x g_ 

2 = cr ■ — — (4) 

|g+xg_| 

can be constructed from the electron polarization vector a and momenta of final state 
particles. If there exists any sizable CP/T violation from EDM of r or in r pair decay 
vertex, in principle the experimental expectation values of these operators are nonzero. For 
EDM of r lepton, d T , the theoretical expectation values of these operators are worked out 



and expressed only as a function of d T \ 13[] . Since at TCF the precision of measurement for 

1CT 3 

these operators are at 10~ 3 level, one expects to probe d T as small as — ~ 1CT 17 e-cm. 

An example is the measurement of d T or d^ in LEP experiment . Expectation value of T* J 



operator is directly related to d T JlOj , 



< T% B >= ^d T C AB diag(-l/6, -1/6, 1/3). (5) 



By the term diag means a diagonal matrix with diagonal elements given above, E cm is the 
energy at cm. frame. The proportional constants Cab depend on the r decay modes, but 
generally this constant is order of one for all the decay models jn|. The decay channels, 



which can be measured in experiments, may be classified as I — I, I — h and h — h classes, 
here I is the lighter leptons, h is charged hadron like tt, p and a\. Very impressively, if the 
initial electron (or both electron and positron) is polarized, one may use the polarization 
asymmetrized distribution. The distribution is defined as the differential cross section dif- 
ference between two different polarizations. With this method, a d T as small as 10~ 19 e-cm 
can be reached at TCF [Tj|, this corresponds to a sensitivity of 10~ 5 of CP/T violation. 
Up to now the best experimental bound on d T is from LEP experimental data, which is 
used to exclude indirectly the d T as large as 10~ 17 e-cm [11], so two order of magnitudes 
improvement on d T measurement can be achieved at TCF. 

Besides the CP/T-odd operator method, several other useful strategy were proposed to 



test these violation in r decay. 1) C. A. Nelson and collaborators [TJ]] investigated systemat- 
ically the feasibility of using the so-called stage-two spin-correlation functions to detect pos- 
sible non-CKM CP violation in the tau-pair production-decay sequence and the correspond- 
ing CP-conjugate sequence. The two-variable energy-correlation distribution I(Ea, Eb,^), 
where \& is the opening angle between the final A and B particles, is essentially a kine- 
matic consequence of the tau-pair spin correlation which depends on the dynamics of Z° 
or 7* — > t~t + amplitude, and of the r~ — > A~Xa and r + — > B + Xb amplitudes. By 
including 9 e and <p e which specify the initial electron beam direction relative to the final- 
state A and B momentum directions in the cm. frame of e~e + system, one obtains the 
so-called beam-referenced stage-two spin-correlation function I(9 e ,(f) e , Ea, Eb,^?)- For the 
7* — ► t~t + vertex, there are four complex helicity amplitudes. Hence, the beam-referenced 
stage-two spin-correlation function constructs four distinct tests for possible CP violation 
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in e e + — > r r + . To illustrate the discovery limit in using the beam-referenced stage-two 



spin-correlation function, Goozovat and Nelson |L5] calculated the ideal statistical errors 
corresponding to the four tests. An advantage of detecting CP violation by use of the stage- 
two spin-correlation function is that the model independence and amplitude significance of 
the results is manifest. It is complementary to the greater dynamical information that can 
be obtained through other approaches, such as from higher-order diagrammatic calculations 
in the multi-Higgs extensions of the SM. 2) Another strategy to test CP violation in the 



two-pion channels of tau decay is due to Y.S. Tsai |16|], the basic ingredient of which is to 
invoke a highly polarized tau-pair. Consider the tau-pair production by electron-positron 
annihilation near threshold. If the initial electron and positron beams are polarized longi- 
tudinally (along the same direction), the tau-pair will be produced mainly in the S'-wave, 
resulting in polarizations of r ± both pointing in the same direction as that of the initial 
beams. Such a polarization is independent of the production angle and the corresponding 
polarization vector supplies us with an important block to form products with the final par- 
ticle momenta. By comparing such polarization-vector-momentum products for a specific 
tau decay channel with those for the corresponding CP-conjugate process, one can perform 
a series of tests for possible CP violation effects in the tau decay. However, it is impossible 
to detect a CP violation in the r — > tiv t decay without violating CPT symmetry. As for 
the two-pion channel, the existence of a complex phase due to the hadronic final-state inter- 
actions, given by the Breit-Wigner formula for the P-wave resonance p, enables detecting 
possible non-CKM violation by measuring asymmetry of (w x qi) ■ q2 without violating 
the CPT symmetry, where (w is the tau polarization vector and q, (i=l,2) are the final 
pion momenta). By limiting the weak interaction to be transmitted only by exchange of 
spin-one and spin-0 particles, one can know that only the S'-wave part of the amplitude for 
the exchange of the extra spin-1 particle make contributions to CP violating observables. 
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A very generic conclusion is that unless two diagrams have different strong interactions 
phases, one cannot observe the existence of weak phase using terms involving w • qi. Tsai 
UlTjl also points out that T violation can not be detected in the pure leptonic decay without 
detecting the polarization of the decay lepton. Because it is impossible to construct T-odd 
operator by the momenta of the initial and final state particles in pure leptonic three body 
decays. This also implies that with CPT symmetry, one can not detect CP violation in 
r decay processes with unpolarized r. On the other hand, however, with polarized initial 
electron and positron beams, one can construct T-odd operators using the momenta and 
polarization vector of r and the decay lepton. Therefore polarization of initial electron and 
positron is very desirable for detecting of CP/T violation at TCF. 3) As for the r — > (37r)v T 
decay, it can proceed either via J p = 1 + resonance a x and the J p = 0~ resonance n'. Choi, 



Hagiwara and Tanabashi |L8| investigated the possibility that the large width-mass ratios 
of these resonances enhance CP-violation effects in the multi-Higgs extensions of the SM. 
To detect possible CP-violation effects, these authors compare the differential decay width 
for the r~ — > 7t + tt~7t~v t with that for the corresponding CP-conjugate decay process. To 
optimize the experimental limit, they suggested considering several CP- violating forward- 
backward asymmetry of differential decay widths, with appropriate real weight functions. 
4) To probe possible CP-violating effects in the tau decay with K~7t~tt + or K~n~K + fi- 



nal states, Kilian, Korner, Schilcher and Wu |TP[ partitioned the final-state phase space 
into several sectors and constructed some asymmetries of the differential decay widths. As 
a result, they showed that T-odd triple momentum correlations are connected to certain 
asymmetries, so their non-vanishing would indicate a possible non-CKM CP violation in 
the exclusive semileptonic r — >three pseudoscalar-meson decays. 

With these knowledge and results obtained in the previous papers in mind, now the 
crucial question, which is also the motivation of this work, is whether for CP/T violation 
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appearing in EDM close to d T ~ 10~ 19 e-cm and CP/T violation effects in r decay as 
10 -3 are possible values theoretically. If for all possible extensions of the SM, which people 
can visualize now, with natural parameter choice, these values are much smaller than the 
theoretically predicted ones, then the effort to search for such small CP/T violation signal 
at TCF would be not much meaningful, at least from the theoretical point of view. In this 
paper we are trying to answer this question by investigating various possible mechanisms 
for generating large EDM of r, CP/T violation in r decay. This paper is organized as the 
following. In section 2 we review the generation of EDM of r lepton in various popular 
beyond standard models and stress on what models can produce possible large EDM of r. 
Following the discussion of EDM, in the section 3 we concentrate on CP/T violation effects 
from t decay in the beyond standard models. The last section is reserved for some further 
discussion, and the conclusion on the possibility of finding CP/T violation at TCF is given. 

2 EDM of r lepton 

EDM of the lepton d\ is a dimension-5 operator. It can only be generated from the loop 
level. Because this operator changes the chirality of the lepton, it must be proportional 
to a fermion mass. In the SM EDM of lepton is generated from three loop diagrams and 



is proportional to lepton mass itself, so di is very small [20|. However generally di can be 
produced from one loop diagrams in beyond standard model. At one loop level, the di can 
be expressed as 

where Mp is some fermion mass, V is a large scale from intermediate states in the loops and 
A denotes other couplings. is a CP/T violation phase. In the following part we assume 
maximal CP/T violation phase, i.e. sin</> ~ 1. From this equation one sees that d\ can 



be at most as large as 1CT 18 — 1CT 19 e-cm if A is between 1.0 — 0.1. Since V is a scale 
around or larger than weak scale, in order to obtain large di, Mp must be a large fermion 
mass such as t quark mass or new heavy fermion masses. For example if M is the r mass 
then d T is smaller than 10~ 20 e-cm which is not detectable at TCF. Same is true for the 
scale V . If V is at TeV scale di is smaller than 10~ 20 e-cm. Although, in principle, d T is 
possibly as large as 10" 19 e-cm , one has to avoid too large EDM of electron d e at the same 
time. Current experimental upper limit on d e is about 10 _26 e-cm. This is a very strong 
constraint especially when one is expecting large d T . So in any beyond standard model, two 
requirements must be satisfied in order to obtain measurable d T . The first one is that the 
model must provide d T at one loop level and d T is not suppressed by a small fermion mass 
term, the fermion mass term should be a top quark mass, supersymmetric partner of bosons 
or other exotic fermion masses. The second one is that the predicted d e associated with 
large d T is below its current experimental bound. These two conditions altogether exclude 
most of beyond standard models which can provide large enough d T observable for TCF. 
We will see from the following discussion that many beyond standard models do not satisfy 
the two requirements. 

Usually EDM of lepton is generated from one loop diagrams in extension of the SM. 
Fig. 1 is a typical one loop diagram for the lepton EDM. The virtual particles are scalar or 
vector boson S and fermion F in the loop. Photon is attached to the charged intermediate 
particles. The di from this diagram is approximately proportional to the fermion mass Mp 
and it is divided by a scale V, which is larger or equal to M F . Besides, there are two more 
couplings at the vertex I — S — F. In a practical model there could be many possible virtual 
bosons and fermions in the loop, but we only consider the dominant contribution here as 
an order of magnitude estimation. The diagram in Fig. 1 is evaluated as 
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where i = e, fj,, r denotes three generation leptons and Q is the electric charge of the virtual 
particles. £ is an order of one factor from the loop integral. Eq. (7) is true up to a factor 
of order one. And there should be a logarithmic dependence on ^ in £, which is slowly 
varying. 

In order to obtain measurable d T and avoid too large d e , one needs a large Mp as 
discussed before and A, A' must be around order of one for r but much small (smaller than 
about 10~ 3 ) for electron. We systematically investigate and review most of the popular 
extensions of the standard model and point out that the following type of models can fulfill 
the requirements. 

Scalar leptoquark mo dels pT]] CP violation effect in r sector for the models are 



recently discussed extensively by some authors JT8|, [22j. It is particularly interesting for 
generating a large d\. These are the models which do not need to introduce additional 
fermion. Because the top quark mass is large, it is possible to generate a large d T through 
coupling of t, top quark and the corresponding leptoquark. d e could be small enough due 
to the coupling of electron, top quark and leptoquark is independent of that for d T . So 
long as there is a relative large hierarchy for the couplings for different generations, the two 
requirements can be satisfied. 

There are five types of scalar leptoquarks which can couple to leptons and quarks. We 
denote them by Si, S 2 , S3, S4 and S 5 . Their quantum numbers under standard gauge group 
transformation are (3, 2, |), (3, 1, — |), (3, 2, |), (3, 1, — |) and (3, 3, — |) respectively. The 
Yukawa coupling terms are therefore given by 

Li = {X\ 3 QLiir 2 E Rj + Ai J t7 ra / Li )Si + h.c. 
U = (Xl J Q^r 2 l c Lj + \2 j UmE c Rj )S 2 + h.c. 
L 3 = \%D Ri l Lj S 3 + h.c. 
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U = \\ J D m EgS 4 + h.c. 
L 5 = {K J QiAir 2 Tl c Lj ) ■ S 5 + h.c. 

Here II and Ql are lepton and quark doublets respectively, U R , D r and E R are singlet quark 
and lepton respectively. Individually only Si and Si contribute to the EDM of lepton. 



£ factor in Eq. (7) is evaluated as £ = |in^f + ^ [p^| . Currently the constraints on 



mass and coupling of leptoquark are relatively weak J23J. For leptoquark coupled only to 
third generation, its lower mass bound is about 45 GeV with order of unit coupling p^. 
This bound is from a leptoquark pair production from LEP experiments. On the other 
hand with the leptoquark mass at weak scale, the coupling is very weakly bounded too. 
In fact the coupling could be as large as order of one. If we take A 33 , A' 33 as 0.5 and the 
mass of leptoquark as 200GeV and assume maximal CP/T violation phase, we estimate 
that d T ~ 2 x 10~ 19 e-cm, while d e is determined by other coupling components, so a small 
d e is not necessary in conflict with a large d T in this model. 

Models with the fourth generation or other exotic lepton The SM with fourth 
generation is another possible model to generate a large d T . The heavy fourth generation 
leptons may play a role of the heavy fermion F in the loop. However it is well known that 
if the fourth generation exists, it must satisfy the constraints from LEP experiments || . 
Here we propose a realistic model for this purpose. 

Besides the fourth generation fermions, we also introduce a right-handed neutrino v R 
and a singlet scalar 77" with one unit electric charge |[25|| . The new interaction terms are 



L = \ijlJiT2ljV~ + KE^urtj- + M R v T R v R + M?v Li v R + h.c. (8) 

where is antisymmetric due to the Fermi statistics. M D is the Dirac neutrino mass 
from standard Higgs vacuum expectation value. In this model three light neutrinos remain 
massless and the fourth neutrino is massive |26|| . The constraints from LEP experiments 
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and other low energy data can be satisfied so long as Mr is at weak scale or up and is 
not much smaller than Mr. In the one loop diagram contribution to d T , t]~ appears as the 
scalar S. The fermion line is two massive neutrinos z/ 4 and vn in the mass basis and they 
are related to each other, 

i>la = cos #z/ 4 — sin 6u H 

Vr = sin Qv± + cos Ovr (9) 

We assume is the lighter neutrino and the dominant contribution is from either vu or 
depending on whether uh is heavier than the mass of 77, M v . d T is evaluated as in (7) with 
M F = M H cos6smd and V ~ M v if M v > M H ; with M F = M U4 cos6smd and V ~ M H 
if M v < M H . Choosing A 34 = A' 3 = 1.0 and M F = 50 GeV, V = 200 GeV, we have the 
numerical result d T ~ 10 -19 e-cm. Also in this model a hierarchy on the coupling A and A' 
for different generation is needed to keep small enough d e , i.e. A 34 >> A i4 and A 3 >> \[. 

Existence of exotic leptons provide another possibility to generate a measurable d T . It 
can be realized in horizontal models |27| . With only three standard leptons, it is impossible 
to obtain large enough d T , because the largest fermion mass in the loop is m T . However, 
with some new heavy leptons this model can provide a large d T . The constraints from low 
energy data can be avoided if one assumes that the horizontal interaction is strong between 
t and the exotic lepton, but it is much weaker in other sectors. Similar result on d T as for 
the case with the fourth generation can be obtained. 

Finally, we should point out that for our purpose it is clear that some new exotic heavy 
leptons are needed in the new physics models, however even though there exists some kind 
of models with some new heavy leptons, they are able to generate d\ only from two loop 



diagrams p8fl , so they may result in interesting d e but not d T . 

Generic MSSM Generic MSSM contains 63 parameters not including the parameters 
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in the non-SUSY SM. Ferminic superpartners of the ordinary bosons can be the heavy 
fermions in the loop diagrams for di. It provides some new sources for CP/T violation. It 
is well known that the electron and neutron can acquire large EDM [[29[ in this model. In 
fact, in order to obey the experimental bounds on d n and d e , some parameters in the model 
are strongly restricted |30| . For d\ generation, it is dominated by photino mediated one loop 
diagram. Both left- and right-handed sleptons also appear in the loop. The contribution 
to d\ from this diagram is proportional to left- and right-handed slepton mixing matrices 
Mlr = (Ai — /itan/3)M;. Ai is the matrix of soft-SUSY-breaking parameters that appears 
in the SUSY Yukawa terms of slepton coupling to Higgs doublet. Here Mi is diagonal mass 
matrix of lepton mass. Usually it is assumed that A\ is diagonal and the diagonal elements 
are not much different for different generation, for example in supergravity inspired model 
Ai is universal for three generation ||, therefore one can get d T /d e ~ m T /m e . Using the 
experimental limit d e < 10~ 26 e-cm, one concludes that d T < 4 x 10~ 23 e-cm []3l |. However 
in the generic MSSM all the elements of A\ are free parameters, so the above constraint is 
not necessarily true. For example if for some unknown reason the 33 component of A\ is 
much larger than other elements, and /i term is much smaller than SUSY breaking scale, 
then d T still can be larger than 10~ 22 e-cm and d e is in the allowed region. In this case d T 
can also be expressed as Eq. (7), but with Mp = m 7 , V = rh^/MLR, A33 = A 33 = e and 
cf) = ar g{M\ n m^) . The loop integral £ was four times the function calculated some years 
ago in dealing with d e in MSSM known as Polchinski-Wise function [52 1. Here m 7 and rh T 
are photino and the third slepton masses respectively. We estimate that d T ~ 10 -19 e-cm 
with m 7 = 100 GeV and V = 200 GeV. 

As for other popular extensions of the SM, we would like to point out here, though they 
have some new sources of CP/T violation, they can not offer a observable d T at TCF. These 



include multi-Higgs doublet model ( including two Higgs doublet model) P, B3|, Left-Right 
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symmetric model [34], mirror fermion model [35] and universal soft breaking SUSY model 



etc. In multi-Higgs doublet model electron [36| and neutron || may obtain a large EDM 



close to current experimental bounds through two loop diagrams, but d T generated in the 
model is quite below the TCF observable value. The reason is that d T is proportional m T , 



but not a large fermion mass. We estimate d T < 4 x 10 e-cm |37| that in this model. For 



Left-Right symmetric model, Nieves, Chang and Pal [38 find that the upper bound for d T 
is 2.4 x 10~ 22 e-cm. It is the right- or left-handed gauge boson in the loop as the role of S 
particle, while right-handed neutrino is the virtual fermion particle in the loop. d T in this 
model is proportional to left- and right-handed gauge boson mixing angle. Though it is 
not suppressed by the small fermion mass ( Mp is a large right-handed neutrino mass), the 



mixing angle is constrained to be smaller than 0.004 [|39| from purely non-leptonic strange 
decays. It leads to about three order of magnitude suppression. In the mirror fermion model, 
standard gauge bosons couple to ordinary lepton and the mirror lepton with a mixing angle. 
It is Z and W bosons in the one loop diagrams, the heavy fermion line is the mirror lepton. 
However the mixing angle in this model is constrained by various experiments flOfl , and 
most stringently by LEP data on Z — > t + t~ The constraint from LEP data on the 

mixing angle is less than about 0.3. The resulting bound is d T < 2.1 x 10 _20 e-cm, which 
is a few times smaller than TCF measurable value. As we have mentioned above in the 
universal soft breaking SUSY model, d T < 4 x 10 _23 e-cm due to the constraint on d e . The 
only alternative situation is discussed above on Generic MSSM in this section. 



3 CP/T violated r decays 



As we have pointed out in the introduction, CP/T violation effects in r decays, if observed, 
must occur at tree level diagrams. That is the interference between the SM r decay processes 
and new tree level processes of r decays, in which CP/T violation phases appear at the 
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interaction vertexes, provides the information of CP/T violation in the r sector. Feynman 
diagrams of these processes can be shown as in the Fig. 2, where f i: fj and fk are light 
fermions. X is a new particle ( scalar or vector boson) which mediates CP/T violating 
interaction. The size of CP/T violation is always proportional to the interference of the tree 
level diagrams. We denote the amplitudes for these diagrams as A 1 for W boson exchange 
diagram, A 2 for other X boson exchange diagrams. The size of CP/T violation in the r 
decay can be characterized by a dimensionless quantity 

\A^+\A 2 \^ [ 1 

Practically physical quantity expectation values which are used to reflect CP/T violation, 
like the expectation values of CP/T-odd operators, difference of a partial decay widths of a 
t~ decay channel and its conjugate r + decay channel, are model dependent and generally 
quite complicate. It needs the detailed information of the new physics model and a lot of 
parameters enter into the expression. This makes it a very much involved work to write 
down these quantities in a specific model beyond the SM. And the exact CP/T violation 
quantity expression written down from a model should be different from the e defined above. 
However as a simple and reasonable estimation, the quantity e in Eq. (11) can be used as 
an indication of how large of CP/T violation may happen at various r decays. Moreover, 
the amplitude A 2 is usually much smaller than A\ because so far all the experimental data 
agree with the SM prediction very well. So A 2 term in the denominator can be neglected. 
Using A\ as the amplitude from W boson exchange and A 2 from the new boson X exchange, 
we estimate its size, 

e ~(4^G F )-M^ (11) 

Here Gf is Fermi constant and A, A' are couplings in A 2 . From Eq. (12) one sees that the 
size of CP/T violation is determined by the parameter - ■ For different models, this 
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parameter is constrained by some other physical processes. So the possible size of CP/T 
violation depends on the parameter region which is restricted in a specific model. 

In Fig. 2 the final state fermions can be a pair of leptons and quarks besides u T . It 
corresponds to pure leptonic and hadronic decays respectively. At the quark level, the 
diagrams with a pair of quarks in the final states denote an inclusive process, it includes 
all possible hadronic channels originated from quark pair hadronization. Some of the useful 
hadronic final states like 2tt, 3tt, Kit, Kirn, KKtt and p, a\ can be used to measure the 
properties of r. However, it is often difficult to make a reliable quantitative prediction 
for CP/T violation in exclusive hadronic decay modes, because of the uncertainty in the 
hadronic matrix elements. On the other hand, for the inclusive cases, one may make a 
more reliable quantitative estimation due to the fact that one has no need to deal with the 
hadronization of quarks in this case. In addition, QCD correction should not change the 
order of the tree level diagram evaluation as the energy scale for r decay processes is around 
lGeV. In this section we only deal with the diagrams containing quark pair inclusively, So 
the CP/T violation size we estimate below is for all the possible hadronic decay channels. In 
the last section we will comment on our results in exclusive processes. Because of the scale 
of r mass, its decay products can only be neutrino, electron, muon and hadrons containing 
only light u,d, s quarks as other heavy quarks are kinematically forbidden. Therefore there 
are not many possibilities for X particle being the candidate for mediating CP/T violation 
in the Fig. 2. In fact all the possible choices are the following: X being leptoquark, charged 
Higgs singlet, doublet and triplet, and double charged singlet. Now we come to discuss these 
different cases separately. 

Scalar leptoquark models At tree level it is obvious that only Si, S 2 and S 5 
contribute to r decays. There are two type of decay processes at quark level, r — > u T ud 
and r — > v T us. The e parameter is determined by A 31 A' 31 * and A 32 A /31 * for these two 
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type of decays respectively in model 1 and 2 in Eq. (8). For model 5 there is CP/T 
violation effect only in the second type process, which is determined by A 32 A' 31 *. A direct 
constraint on these parameters can be obtained through comparing the theoretical value 
T th (r -> txv t ) = (2.480 ± 0.025) x 10~ 13 GeV and the measurement value of T exp (r -> 
ttu t ) = (2.605 ± 0.093) x 10~ 13 GeV |42|]. Assuming that real and imaginary part of the 



coupling AA'* are approximately equal, one has from r — > ixv T [FJ 
|Jm(A 31 A' 31 *)| |i2e(A 31 A /31 *^' 



< 3 x 10~ b GeV (12) 



at 2a level for model one and two. And from r — > Kv T a similar result can be obtained for 
all the three models. Using the theoretical value T th (r —> Kv T ) = (0.164 ± 0.036) x 10~ 13 
GeV PL H and the measurement value r exp (r -> Kv T ) = (0.149 ± 0.051) x HT 13 GeV for 



the r — > Kv T decay width we obtain 

|/m(A 32 A' 31 *)| |ii!e(A 32 A /31 *) 



M\ M x 



< 7 x l^GeV (13) 



at 2cr level. This constraint is less stringent due to the large uncertainties in T exp (r — > Ku T ). 
With these constraints, one estimates the upper bound of e value for the two type of processes 
as 

e ( r - ^ Vr ud) ~ (4v / 2G F )' l Jm(A . 3 . 1 2 A/3U) < 4 x 10- 2 (14) 

M x 

and 

e(r" -> i/ T iis) ~ (4v / 2G F )- 1 sin 6 C ^^ - < 2 x 10~ 2 (15) 

where #c is Cabibbo angle. e(r~ — ► z/ r us) is proportional to sin 6^ and is smaller than 
e(r — > u T ud) because this process is Cabibbo suppressed, even though the coupling is less 
constrained than that of Cabibbo unsuppressed process. From this estimation we expect 
CP/T violation in these models could be large enough for TCF or in the other words 
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TCF data can put stronger direct restriction on the parameters of the model. However, 
if one assumes that all the couplings A and A' are at the same size irrespective of the 
generation indexes, then much more stringent bounds exist. These bounds are obtained from 
experimental bounds of BriK^ — > fie), Br(ir —>■ ei/ e (j)), Br(TT — > fiv^)) and T(fiTi — > 
eTi)/T(fiTi capture) ||18|| . They are generally about five order of magnitude smaller than 
the direct bounds. Therefore the size of CP/T violation is e < 4 x 10 -7 which is far below 
the capability of TCF. 

Multi-Higgs doublet models (MHD) With the natural suppression of flavor 
changing neutral current, it is necessary to have more than two Higgs doublets, so that 
there are at least two physical charged Higgs particles. CP/T violation may generally 
happen through the mixing of these charged Higgs particles. We consider a multi-Higgs 
doublet model, say, n Higgs doublets. In this model there are 2(n-l) charged and (2n-l) 
neutral physical scalars. Since only the Yukawa interactions of the charged scalars with 



fermions are relevant for our purpose. Following Grossman |43J] we write down the Yukawa 
interactions in fermion mass eigenstates as 



'MHD 



2V2G F ^ =2 [X i (U L VM D D R ) + Y^UrMuVDl + Zi(J L M E E R )]Ht + h.c. (16) 



Here M V) M D and M E denote the diagonal mass matrices of up-type quarks, down type 
quarks and charged leptons respectively. V is KM matrix. X, Y and Z are complex 
couplings which arise from the mixing of the charged scalars and CP/T violation in r 
decay processes is due to these couplings. How large is the e for various r decay channels 
depends on the values of these parameters. More precisely, in the pure leptonic decays 
the size of CP/T violation is determined by Im(ZiZ*) with i ^ j and in hadronic decays 
it is determined by Im(XiZ*) and Im{YiZ*). The three combinations of parameters are 
constrained by various experiments jPJ. The strongest constraint on Z is from e — // 



universality in r decay, which gives \Z\ < 1.93M#GeV 1 for Higgs mass Mh around lOOGeV. 
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Im(XZ*) is bounded from above from the measurement of the branching ratio Br(B — > 
Xtv t ), Im(XZ*) < \XZ\ < 0.23M^GeV~ 2 if M H < 440 GeV. Finally a upper bound is 
given as Im(YZ*) < \YZ\ < 110 from the experimental data of the process K + — > 7r + z/z/This 
bound is obtained for t quark mass at 140 GeV |44| and Mh = 45GeV, however for a different 
Mhi say lOOGeV, this bound is expected not to change much. With these bounds we can 
estimate CP/T violation size of r leptonic and hadronic decays. For the leptonic decay 
t — > fivv, we have the quantity 

ItamVV . m, = l^ Im{zz . } < 2 x (17) 
2 M 2 H m T 2M 2 H y ' ~ y ' 

Here the additional factor — comes from the interference of left- and right-handed muon 

lines in the final states. So we expect that CP/T violation effect in the process r — > evv is 

suppressed by a factor rrie/m^ and is negligible. For the hadronic decay r — > -urfz/ we have 

1 m d m d 



2 M| 



-Im(XZ*) < 3 x 10" 4 , (18) 



With the current d quark mass = 7 MeV and the dynamical d quark mass = 300 
MeV . For hadronic decay r — > usv similar result is obtained 

e - 2^T /m(xz } - L5 x 10 (19) 

Here we use current and dynamical s quark masses as 150 MeV and 400 MeV respectively. 
In summary, in multi-Higgs doublet model CP/T violation effect is possibly as large as order 
of 10 -3 for exclusive hadronic decays and It could be even close to 10~ 2 in pure leptonic 
decay to \i and neutrinos. 

Other extensions of the SM for pure leptonic decays Besides leptoquark 
and Higgs doublet, there are three other kind of scalars which can couple to leptons. We 
denote I as a lepton doublet and E as a singlet lepton. Two / can combine to a charged 
singlet or a triplet. Two E can combine to a double charged singlet. Corresponding to 
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these three cases one can introduce a charged singlet scalar h~, triplet scalar A and double 
charged scalar K . However K only induce a lepton family- number- violating process 
t — > 3/. There is no diagram corresponding SM contribution, so there is no CP/T violation 
mediated by this particle. Also the branching ratio (< 10~ 5 ) for this decay is much smaller 
than TCF reachable CP/T violation precision 10~ 3 . In principle if there exists more than 
one h or A, CP/T violation can be induced by the interference of the W exchange diagram 
and h or A exchange diagram in the process r — > Wu with I = e, fi. Now let us discuss these 
two possibilities in details. We can write down the new interaction terms which couple the 
new scalar particles to leptons as the following 

Lh = \fifiCiT2lih + h.c. (20) 

La = -gijfiCiTzrljK + h.c, (21) 

where C is the Dirac charge conjugation matrix and is antisymmetric, is symmetric 
due to Fermi statistics, e parameter for these singlet and triplet models are given by, 

e h * (4V2G F )^ Im{ ^ l) (22) 



in singlet model and 



-i Im{g r igl T ) 



e A * {AV2G F rt^IlL (23 ) 
1V1 A 



in triplet model respectively. 

For the singlet model we assume that f efl is considerably smaller than f T i, so that one 
does not need to readjust the Fermi constant Gp- This assumption is also consistent with 



the constraint set by universality between j3 and /x decay |25], ^5[. The parameter Im ^-± lT) 



is constrained only by the measurement of r leptonic decays. At 2a level (which is about 
2 ~ 3% precision) we estimate approximately lm ^p lT ^ < 1CT 6 GeV~ 2 f46|. It implies that 

h 



-Jm{f Tl ft T ) 



t h ~ (W2GV)" 1 Vrl < lA x 10_2 (24) 
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with Mh = 100 GeV. Therefore in this model there is a possibility that CP/T violation 
effect may show up with a size reachable at TCF in pure leptonic decay channels. 

For the triplet model the direct constraint is also from the measurement of pure leptonic 
decays. The same result is obtained as that in the sing let model , i.e. < 10~ 6 

h 

GeV -2 . As the result of this constraint one has 



eh ^ 



(AV2G F )- l Im [ 9 ;f r) < 1.4 x 10- 2 (25) 



Mi 



with Ma = 100 GeV. However, in this model the new interactions will induce lepton family 
number violating decay r — > 3/ and \x —>■ 3e through exchange of the double charged scalar 
particle A . Without seeing any signal, one obtains some approximate bounds on the 



coupling constants as the following [47 



eel < 5 x 10~ 12 (26) 



Ml 



and 

\9t19u\ 



Mi 



< 10~ 8 (27) 



for Ma = 100 GeV. If one assumes that all the couplings <7y are at the same order of 
magnitude, then these bounds will restrict the CP/T violation size far below the ability 
of TCF. Again we see some hierarchies on the couplings are needed for this model to give 
rise observable CP/T violation effects. Additionally in the triplet model one has to avoid 



the restriction from neutrino mass generation |48|. If neutrino develops a mass at tree 
level, either the couplings or the vacuum expectation value of the neutral component of the 
triplet A are extremely small. The natural way to deal with this problem is to impose 
some symmetry on this model. An example is to introduce a discrete symmetry: 

il; E -+iE A -> -A (28) 
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With this symmetry, A will never develop a nonzero vacuum expectation value, therefore 
the couplings are not constrained by the neutrino mass generation. 

4 Discussion and conclusion 

In this work we systematically investigate the possibility of finding CP/T violation in the 
r sector with TCF. The origin of CP/T violation is from the extensions of the SM. We 
discuss most of the popular beyond the SM and present the models which may give rise 
large CP/T violation in r sector through either EDM or decay of r lepton. Before making 
our conclusion, some interesting points should be further discussed or emphasized. (1) 
Polarization of initial electron and/or positron is very desired for our purpose. First with 
polarization the precision of measurement of EDM will be increased by about two order of 
magnitude, as lCT 19 e-cm, which is used through this work. Without polarization, from our 
above discussion one sees that we have no hope to expect a detectable EDM of r at TCF. 
Secondly in some decay channels without final state interaction, like pure leptonic decays 
and two body decays nv T etc., polarization is needed to search for CP/T violation occurring 
at t decay vertex. With unpolarized electron and positron beams the CP/T violation could 
only be detected using channels with final state interaction phase, like 2irv T etc. (2) For the 
hadronic decay we only consider inclusive processes. The advantage of inclusive process is 
that one does not need not to consider the hadronization of quarks, which may bring in large 
uncertainties in the estimation. And the event number in inclusive process is larger than 
that in certain exclusive processes. However we should mention that for certain exclusive 
decays the CP/T violation parameter e can be larger than that in inclusive decay. One 
example is from the multi-Higgs double model. We estimate that e < 3 x 1CT 4 for the decay 
r — > udu T . Here we may also consider the exclusive decay r — > ?>-nu T contributed by a\ and 
7r' resonances. Compared to inclusive decay, the e parameter is larger by a factor of (using 
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current algebra relation) 

< o\uLdji\n' > m n i 



100. (29) 



< o\u L nf d L \ir' > m u + m d 
So e < 3 x 10~ 2 is obtained. However on the other hand the event number decreases by a 
factor of 

fl Br(r — > 7TU T ) 9 

±2L l II ~ iq- 2 (30) 

f n Br{r — > hadron + v T ) 

Here = 5 x 10 3 GeV is used. Therefore statistical error increases by about 10 times. In 
the other words the measurement precision at TCF for this channel is about 10~ 2 . As the 
result, at 2a level e ~ 3 x 10~ 2 is observable. This estimation agrees with the exact result of 



reference [|18|]. (3) Obviously the numerical result we obtained above is quite crude. More 
accurate estimation is necessary in the future. For instance through this paper we assume 
that EDM as large as 10 _19 e-cm and e as large as 10~ 3 can be observed. This of course is a 
rough estimation. To be more precise, Monte Carlo simulation is needed, which will tell us 
more confidently how large CP/T violation is able to be observed at TCF. Especially the 
Monte Carlo simulation on EDM of r will give us a quite clear result , because in this case 
the d T is the only parameter we should take care. All the model dependence is included in it. 
Recently a group of people analyzed the data from BEPC experiments to set bound on the 
T- violating effect for r system ||49[| . Following the suggestion by T.D. Lee, they considered 
the pure leptonic r decays to e fi^ plus neutrinos in the final states. The T-violating 
amplitude 

A=<p e - {p x X p 2 ) >average (31) 

is measured, where p e is the unit momentum vector of the initial electron beam, pi and p 2 
are the unit momenta of the final state electron and muon respectively. Totally 251 events 
are analyzed and it results in 

A = -0.097 ± 0.039 ± 0.135 (32) 
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This result agrees with no T-violation as expected from our previous discussion on pure 
leptonic r decays. (4) In order to generate detectable large CP/T violation effects, we know 
from our investigation that there must exist new physics and the new physics scale is not far 
above the weak scale. Therefore if there is a observable CP/T violation effect in r sector at 
TCF, the associated new physics phenomena should be observed at high energy experiments, 
like LHC and LEP II experiments. It is interesting to see if the new particles predicted by 
the various models we have discussed in this paper are indeed detectable in these high energy 
experiments. (5) Precise measurement of the pure leptonic decay is another way to test the 
new physics responsible for CP/T violation. Since if there is CP/T violation effect at level 
of 10~ 3 , the r leptonic decay width must deviate from the SM prediction at the same level. 
So we expect to observe the deviation by measuring the branching ratio of the pure leptonic 
decay. However it is not true vice versa, since a deviation of leptonic branching ratio from 
that of the SM does not necessarily indicate CP/T violation. 

Finally we come to our conclusion. There exists the possibility that CP/T violation in 
r sector is large enough to be discovered at TCF, although for this large violation effect 
some specific new physics phenomena beyond the SM are needed and the parameter spaces 
of the models are strongly restricted. 
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Figure Caption 

Fig.l One loop diagram for lepton EDM generation, where F is heavy fermion and S is 
the new boson. Photon line is attached to charged particles in the loop. 

Fig. 2 The diagrams for r decay, (a) is the contribution from the SM and (b) is the 
contribution from new boson exchange. 
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